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Abstract 

In this paper, we present an implementation of the harmonic polylog- 
arithm of Remiddi and Vermaseren [1] for Mathematica. It contains an 
implementation of the product algebra, the derivative properties, series 
expansion and numerical evaluation. The analytic continuation has been 
treated carefully, allowing the user to keep the control over the definition 
of the sign of the imaginary parts. Many options enables the user to adapt 
the behavior of the package to his specific problem. 
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1 Introduction 



The harmonic polylogarithms (HPL's) [1] are a generaUzation of the usual poly- 
logarithms [2] and of the Nielsen polylogarithms [3]. The HPL's appear in many 
calculations in high energy physics. They are found in three-loop deep inelastic 
splitting and coefficient functions [4-7], in two-loop massive vertex form fac- 
tors [8-15], in two- loop Bhabha scattering [16-21] and in multi-loop three-point 
and four-point functions [22-32]. The HPL's also show up the expansion of 
hypergeometric functions around their parameters [33-36] and in more formal 
developments [37]. 

The HPL's have already been implemented for the algebraic manipulation 
language FORM [38] . This implementation is described in [39] and can be found 
at [40]. A FORTRAN code for the numerical evaluation of the HPL's is also 
available [41]. A numerical implementation is also available for GiNaC [42], 
described in [43]. The aim of this work is to provide an implementation of the 
HPL's, including numerical evaluation. A Mathematica implementation of the 
HPL's will in particular allow the usage of these functions in connection with 
many of the multi-purpose features of Mathematica. Concerning speed, it is 
not expected that our Mathematica implementation outperforms the FORM 
implementation [39]. 

The paper is articulated as follows. In section El we review the definition of 
the HPL's. Their analytic properties are described in section 13 (derivative) , 01 
(product algebra) and |S1 (singular behavior). In section Owe treat the different 
identities relating different HPL's. Section [71 presents the series expansion of 
the HPL's. The HPL's of related arguments are treated in section |S1 Section 
El treats the analytic continuation of the HPL's. Values of HPL of argument 

1 are related to the multiple zeta values (MZV), which we discuss in section 
Cni In section [TTl we describe our implementation and provide examples of its 
usage. Finally we conclude in sectional Some tables and calculational details 
have been relegated to the appendix. In appendix ^ we list some identities 
for the MZV. In appendix IbI we prove the relation between the HPL's of unity 
argument and the MZV's. Finally in appendix ICl we list the representation of 
some HPL's in terms of more common functions. 

2 Definition 

The harmonic polylogarithms (HPL) H (ai, . . . , a/t; a;) are functions of one vari- 
able X labeled by a vector a — (oi, . . . , a^). The dimension k of the vector a is 
called the weight of the HPL. We define the functions 
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/-i(^) = t]- (1) 
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The HPL's are defined recursively through integration of these functions. For 
weight one we have 



H{l-x) 


= 1 h{t)dt = J 


^dt = -log(l-x) 









H{0;x) 


= log(a;) 

X 


X 


Hi-hx) 


= y/_i(i)di= 


1 Y^dt = log(l + a;), 
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and for higher weights 



X 

H{a,ai^...X,x) = j fa{t)H{ai^...x,t)dt, (3) 



where we used the notations 

"'i = i,...,i, and ai,...,fe = ai, . . . , a*;. 

n 

A useful notation introduced in [1] for harmonic polylogs with non-zero right- 
most index is given by dropping the zeros in the vector a, adding 1 to the 
absolute value of the next right non-zero index for each dropped 0. This gives 
for example (3, —2) for (0, 0, 1, 0, —1). We can extend this notation to all index 
vectors by allowing zeros to take place in the extreme right of the new index 
vector. This gives for example (3, —2, 0, 0) for (0, 0, 1, 0,-1, 0, 0). We will enclose 
index vectors in this notation in curly brackets and refer to it as the "m"- 
notation, as opposed to the "a" -notation. Some formulae or transformations 
are easier expressed in the one or the other notation, therefore we keep both 
notations in parallel. 



3 Derivatives 

The formula for the derivative of the HPL's follows directly from their definition 

-^H{a, ai,...,fe; x) = fa{x)H{ai^...^k]x). (4) 
In the "m" -notation, the derivation reads for n positive 

-^H {{n,mi^,„^k};x) = l,mi,...,fc};a;), n>l 

-^H {{-n,mi^...^k};x) = ^-^({-(n- l),mi,...,fe};a;) , n > 1 

-^i? ({l,mi,...,fc};a;) = -^—H{{mi^...^k};x) 
ax 1 — X 

{{-l,mi^...^k};x) = Y^i? ({mi,...,fe};x) 

^H{rO};x) = x) (5) 

4 Product algebra 

The product of two HPL of weights wi and W2 can be expressed as a linear 
combination of HPL's of weight w = wi + W2- The formula in the "a"-notation 
for two HPL's with index vectors p and q is given by 

H{pi,...,p^^;x)H{qi,...,qy,^;x) = H{p;x)H{ci;x) = ^ H{r;x) 

repttlq 

where p W q is the set of all arrangements of the elements of p and q such that 
the internal order of the elements of p and q is kept. For example one has for 

p = (o, b) and q = (y, z) 

H{a,b;x)H[y, z;x) = H{a,b,y, z;x) + H{a,y,b, z;x) 

+H{a, y, z, b; x) + H{y, a, b, z; x) 
+H{y,a,z,b;x) + H{y,z,a,b;x) (6) 
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5 Extraction of the singular behavior 



The HPL's can have divergences in a; = and x = 1. The divergent part can be 
extracted with the help of the above product rules. Divergences in appear as 
logarithmic divergence with behavior ~ log" (a;) if there are n O's at the right end 
of the index vector. In this case one can make the divergent behavior explicit 
by writing 

-ff(ai,...,fc,0;a;) + H{ai^,,,^k-i,0,ak;x) H h i?(0, ai, x). (7) 

which one can solve for _ff(ai^...^fc, 0; x). Recalling that H{0;x) = log(a;), the 
divergent log stands now explicitely. If there are more zeros in the right end of 
the index vector, one has to use this method recursively, until all end zeros have 
been exchanged against log's. 

Similarly, when there are n I's at the very left of the index vector there 
appears a logarithmic divergence ~ log"(l — x). This can be extracted in the 
same way as above. One uses the product rule 

H{l;x)H{ai^,,,^k;x) = 

-ff(l,ai,...,fc;a;) + H{ai,l,a2,...,k;x) h i?(ai,...,fc, 1; x). (8) 

Here we solve for _fr(l, ai^...^^; x). Again we can apply the method recursively 
to make the singular behavior by x = 1 explicit. 

This method allows to express all HPL with left I's and/or right O's in terms 
of products of H{l;x), H{0;x) and HPL's without left I's and right O's. 

Using the same method, it is also possible to extract simultaneously the 
H{—1; x) and H{0; x) coming from the —1 on the left and on the right of the 
index vector. The factors 1; x) will correspond to divergences in —1 which 
are relevant for the analytical continuation to the interval (— oo, —1), see section 
El One cannot however extract simultaneously the H{—l;x) and H{l;x) by 

iJ(l,-l;x) = iJ(-l;x)ff(l;x) - iJ(-l,l;x). 

We will make extensive use of these relations in the next sections. 

6 Minimal set 

The procedure for extracting the divergent parts of an HPL described in section 
|S1 allows to express many HPL's in terms of HPL's without divergences at 
and 1 (so called "irreducible") and products of HPL's of smaller weight. The 
product rules described in section 01 also provide relations between HPL's of a 
given weight and HPL's of lower weight. One can combine all these relations 
to get a minimal set of HPL's for a given weight from which one can construct 
all other HPL's of this given weight, up to products of HPL's of lower weight. 
Tabled shows the number of HPL's, irreducible HPL's and the dimension of 
the minimal set as a function of the weight. 

Only the number of elements in the minimal set is fixed, there is a freedom left 
for the choice of which elements are to be taken as independent. Our choice 
was first to exclude all HPL's whose divergent behavior can be extracted along 
the lines of section jSl from the minimal set. For the remaining ones, we ordered 
the index vectors with the following procedure. One adds one to all indices of 
the index vector (in the "a" -notation), the result is to be interpreted as the 
expansion in basis 3 of a number. This number describs the index vector on a 
unique way. We used this numbering to sort the irreducible HPL's and choose 
to express the last as a function of the first. 
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Weight 


Full basis 


Irreducible set 


Minimal set 


1 


3 


3 


3 


2 


9 


4 


3 


3 


27 


12 


8 


4 


81 


36 


18 


5 


243 


108 


48 


6 


729 


324 


116 


7 


2187 


972 


312 


8 


6561 


2916 


810 



Table 1: Dimension of the different basis 

7 Series expansion 

The series expansion of the HPL's can be defined recursively. Let us call 
Zi{ai....^k) the coefficients of the expansion of H (ai....^fc; x). 



(9) 



We assume that the index vector has no trailing 0, as these lead to log(x) 
divergences. These divergences have to be made explicit by the procedure of 
section 121 Now we use the definition of the HPL's 



nX nX 

/ dx'fa{x')H{ai,,„x,x')^y2Zdai,...,k) dx' fa{x')x' 
Jo Jo 



For the three different possibilities 1, 0, —1 for a we get 



dx' II 

X 



l-x' 
" dx' 

X 

dx' 



pX '-^ -1 







dx' ,1 
7^ 



X 

i 



1 + 



/ dx'{xr Y.^-xy = (-1)^+1 E 



J 



(10) 



The recursion relation for the Z^'s is found by interchanging the order of the 
summations over i and j 



Zj (-1, ai,...,fe) 



1 ^ 

-y^^i-i(ai,...,fc) 



i=2 



— -^j"(ai,...,/c) 



E(-ir+'^.-i(ai,....fc) 



(11) 



For the "m" -notation, we can use the same notation for the coefficient of the 
series expansion, as no confusion is possible. 

oo 

Here again we assume that all log(a;) divergences have been made explicit by 
the procedure of section O so that the vector mi^,,,^k has no trailing 0. The 
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recursion relations for n positive reads 

1 ^ 

Zj{n,mi^,,,^k) = — ^Zi_i(mi,...^fc) 

Z,i~n,m,_k) = ■ti^^(-iy+iZ,_i(mi,...,fe). (12) 
The recursion is settled by the expansions of H{{n}; x) and H{{~n}] x) 

OO j 

H{{^nYx) = -iz„(_:,) = ^^t-^, (13) 



which give 



sgn(n)'+i . ^ „ 

0, i < . 



8 HPL's of related arguments 

In this section, we present the identities between HPL's of different (but related) 
arguments. All the identities are obtained through change of variable in the 
definition of the HPL's. The easiest case is the change of variables 

X — —t. 

This transformation is only allowed if the right-most index is not zero, as in this 
case one would get negative argument in a logarithm^. 

i/({mi,...,TOfc},-a;) = {-\f H {{-m^, . . . ,-vn^}\x) (14) 

This identity holds in the "m" -notation. For index vectors in the "a" -notation, 
the exponent k is not the length of the vector, but the number of 1 and —1 in 
the index vector. 

We now consider the change of variable 

X x^ . 

Since we cannot express 1 + x"^ as sum of our basis functions /i,/o and 

we will exclude index vectors with negative indices for our considerations. The 

identities for the weight 1 are 

H{0;x^) = \og{x^) = 2H{0;x), 

H{l-x^) = \og{l- x^) ^ H{l;x) - H{-l]x). (15) 

There we chose x to be the positive root of x^ . Furthermore we see in the second 
equation that the identity only holds for x smaller than 1. For higher weights 
we use recursively the relations 



/"^^ dx' 

H{Q,m2...,k]x^) = / —H{m2 k]x') 

Jo 

r^^ da;' 

H{l,m2...,k\x^) = / -H{m2,...,k;x') 

. n l-x' 



fdt'(^-^)i7(^,...,;t'% (16) 



'^See section |n] for the analytical continuation in this case 
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where H{m2....,k', i'^) is expressed as HPL's of argument t' , which is known in a 
recursive approach. 

The next transformation we consider is 

X I — X. 

Since 1/(2 — a;) (the transform of 1/(1 + x) — can not be expressed as hnear 
combinations of the /i's, we will only consider index vectors without negative 
index. We first have the identities for weight 1 

H{0;l^x) = log(l-2;) = -i/(l;a;), 

H{l;l~x) = -log(a;) = -iJ(0;a;). (17) 

These identities only hold for x between and 1. Here again, we process recur- 
sively in the depth. We use the fact that one can express HPL's with I's on the 
left of the index vector as product of H(l;x) (whose transformation we know 
from above) and HPL's without left I's and treat only the case of a as the left 
index. For this we use the formula 



f^^"^ Ax' 

I —H{m2,...,k;x') 
Jo ^ 

= / -n-H{m2,...,k]x)- I —H[m2,...y,x) 

r dt' 

= i/(0,7ni,....fc;l)- / -H{m2,...,k;l-t'), 

Jo -L — I 

(18) 

where one has to insert for H{m2,...^k] 1 — t') the expansion in terms of HPL's of 
argument t' known from the recursion. The HPL's (with positive m's) evaluated 
for argument 1 are related to the multiple zeta values (MZV). This aspect is 
treated in section EH 

We consider now the transformation^ 

1 1 

x = — hie, t= ze, x > \. 

t X 

The identities for weight 1 read 

H{Q;x) = -H{0;t), 
H{l;x) = H{l;t) + H{0;t) + iTT, 
H{-l-x) = Hi-l;t) ^ H{0;t). (19) 

For higher weight we proceed by induction. We will assume that the leftmost 
element in the vector index is or —1, since we can extract the one's with the 
procedure described in sectional 

X 

H{a,a2,...,k]x) = f dx' fa{x')H{a2.,...,k]t) 

dx' - j dx' \ fa{x')H{a2,...,k;x') 

X / 

1 

H{a, a2_k; 1)+J ^fa (^^) H (^a2,...,fe; ^) (20) 



^The convention for the sign of the imaginary part of x is different from [1]. 
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We now distinguish the two cases a = and a = — 1, 

1 



Inserting for 



its representation in terms of HPL's of argument t' (which is known in a recursive 
approach) one gets terms of the kind 

1 

J dt'h{t')H{b2,...,k; t') = H{b, b2,...,k; 1) - H{b, b2,...,k;t). (22) 
t 

The term H{b, 62....,fc; 1) is a finite constant, as b is either —1 or and not one. 
The next transformation is 

x = t/{t-l). 

Since log((2t - l)/(i - 1)) (the transform of log(l + x)) and log{t/{t - 1)) (the 
transform of log(a;)) can not be expressed in terms of fifl,-i{t), we restrict the 
following consideration to HPL without trailing O's and without negative index. 
We have only one identity of weight 1, 

H{l;x) = -H{l;t) , 

which only holds fov x,t < 1. For the higher weights we use the recursion 
relations 



H{0,m2...,k;x) = / —H{m2,...,k;x') 
Jo 



'dt' (--^]H(m2,...,k; ^' 



H{l,m2...,k;x) = / —H{m2,...,k]x') 
Jo x' 



where the HPL in the integrand is replaced by its expansion in terms of HPL's 
of t, as known from the recursion. 

The last transformation we consider is 

1-t 
X = . 

l + t 

Here wc can also consider negative index in the index vector. The identities for 

weight 1 are 

H{Q;x) = -H{l;t)-H{~l;t), 
H{l;x) = -H{0;t)~H{-l;l) + H{-l;t), 
H{-l;x) = -Hi-l;t) + Hi-l;l). (24) 
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For the identities of higher weight, one uses again the fact that HPL's with I's 
on the left of the vector index can be reduced. The recursion relations for the 
two other cases are 



/ —H[m2 fe;x) 

= / —H{m2,...x,xj- —H{m2,...,k]x) 

JQ X Jx X 

= H{Q,m2,...,k\l) 



H{-l,m2...,k;x) = 



/ —H{m2,...,k]x) 
Jo ^ 

f* d<' 

HiO,m2^...,;l)^l^H 



+ t' 



1-t' 

7719 

-•• ■ ' 1 + t' 



where we insert for 
the identity of lower weight known from the recursion. 



H f m2,...,fc 



(25) 



(26) 



9 Analytical continuation 

The transformations of the preceding section were valid for restricted intervals. 
In this section we consider the analytical continuation to the remaining real 
axis. We first consider the continuation to the interval (—1,0). We define 

X = —t + Sei , 

where < t < 1 and e is infinitesimally small and positive. 6 is either +1 or 
— 1. The continuation for the HPL of weight 1 are 

H{l;x) = H{-l;t), 
H{0;x) = H{0:t) + Si7r, 
H{-l-x) = H{l;t). (27) 

For higher weight one extracts the factors H{0; x) with the method of section 
Eland substitutes the above result for it. For the remainder, one can use the 
formula fTl| of the preceding section. 

The next interval to consider is the interval (l,oo). There we use the trans- 
formation 

1 

X ^ t = — , 

X 

which transforms the real axis greater 1 into the interval (0, 1). We consider 

x = — \- iSe , i6e , a; > 1. 

t X 

This transformation has already been treated in the preceding section for S — I. 
We can use the same procedure to keep track of the sign of the infinitesimal 
imaginary part by modifying the transformation of H{l]x) to 

H{1; x) = H{l-t) + iJ(0; t) + iSi: . 
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The last interval is (— cxd, —1). There one also uses the transformation x —>■ —x. 
Again we consider 

x = -t + 5ei, t^-x-Sei. (28) 
The identities of weight 1 are 

H{l;x) = -H{-l;t) = -H(-l:]:]+H(o;l 



H{0-x) = H{0;t) + Sin = -H[0;-] +6iTT 



1 

t 



H{-l;x + Sei) = -H{l;t - Sei) = 

For higher weight one extracts both the H{0;x) and the H{—l;x) with the 
method of section [S] and substitutes the above expressions. The remainder is 
transformed with 1)14(1 . Since t is larger than 1, one should use the analytical 
continuation described above to bring the argument of the HPL into the interval 
(0,1). 

10 Values at unity and Multiple Zeta Values 

Some transformations of the preceding section, there appear HPL's of argument 
1. These are related for positive m's to the Multiple Zeta Values (MZV) and 
for general m's to colored MZV. The relation can be found through induction"^ 
and reads 

H{{mi_k};l) = ^(mi,...,fc)C(mi,...,fc), k>l (30) 
H{{m};l) = Cim), m>0 (31) 

H{{-m};l) = (1 - 2i-'")C(m) m>0. (32) 



The MZV C are defined as 



oo 21—1 "ifc-i^l k 



am.,...,m,)^j:E--- E n^^l3^- (33) 



■ \mj\ 

il=l 1 1 3 = 1 

and are described in the literature, for example in [44,45]. The vector m is 
obtained from the vector m with 

m = (mi, sgn(mi)m2, . . . , sgn(mi_i)mi, . . . , sgn(mfc_i)mfe) (34) 

The factor N{mi^,,,^k) is given by 

N{mi_k) = (-l)#("-<o) (35) 

The MZV's also form an algebra. Due to this fact, they can be expressed in 
terms of a few mathematical constants like powers of tt, ^-functions and polylogs 
at specified values. We list some of the identities of [44,45] in Appendix 1X1 For 
the implementation of the HPL at unity, we translated the tables of the FORM 
package harmpol.h [1] and their expansions for weight 7 and 8 htable7.prc 
and htableS.prc for Mathematica. 

In these tables, there appear some constants that are not expressible through 
known constants like tt, C{n), log(2), or Li„(l/2). Using the different rela- 
tions between the different MZV, one can reduce the number of independent 

^see Appendix IbI 
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constants. Which constants are kept is a matter of choice. In the tables we 
translated, the choice of the independent constants is 



S6 



S7a 



ssb 



il=l 1 i2=l 



C(-5, -1) + C(6) ~ 0.98744142640329971377 (36) 



21=1 12 = 1 13 = 1 

C(-5, 1, 1) + C(-6, 1) + C(-5, 2) + C(-7) (37) 
-0.95296007575629860341 (38) 

1 ^ (-1)*^ ^ (-1)^3 



il— 1 ^ 22 — 1 ^3— 1 

C(7) + C(5, 2) + C(-6, -1) + C(5, -1, -1) (39) 
1.02912126296432453422 

CO ^ ii ^ 



ssa = 5({5,3},oo) = ^ ^ ^ - 



2" ^ * 1^ 

1 12 = 1 2 



C(8) + C(5, 3) ~ 1.0417850291827918834 

(-1)^1 (-1)'' 



5({-7,-i},oo)=5:i^5: 



11=1 12=1 



= C(8) + C(-7, -1) =i 0.99644774839783766600 (40) 
S8c = S'({-5,-l,-l,-l},cx3) 

il = l 12 = 1 23=1 «3 = 1 

= C(8) + C(-7, -1) + C(-5, -3) + C(6, 2) + C(-5, -1, 2) 

+C(-5, 2,-1) + C(6, -1, -1) + C(-5, -1, -1, -1) (41) 
~ 0.98396667382173367094 

°° / iWi »i / -|^^2 *2 , »3 . 

- s((-5,-i,i,i),»)^E^E4^E^Et 

i 1 — 1 ^2 — 1 ^3 — 1 ■^3 — 1 

= C(8) + C(-5, -3) + C(6, 2) + C(7, 1) + C(-5, -2, 1) + C(-5, -1, 2) 
+C(6, 1, 1) + C(-5, -1, 1, 1) 0.99996261346268344768 (42) 

For index vectors with left I's, there appear divergences of the form 

oo 

^- = 5(l,(^) = i/(l;l). 
1=1 

These divergences are well defined^ and can cancel during a calculation. Right 
O's can be extracted with the method of section[S| All factors -ff (0; 1) = log(l) = 
vanish, even multiplied with iJ(l; 1) = — log(O), since 

log(a;) log"(l - x) ^ 0, a; ^ 0, n > 0. 



11 Implementation 

In this section we present our implementation of the HPL's. The package can 
be found at [46]. Installation instructions can be found there. After instalation, 
the package can be called with 

"see [1]. 
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«HPL' . 

This call should be done at the beginning of the Mathematica session. 
11.1 New functions 

In the package HPL we defined the following new functions 

• HPL [m,x] is the harmonic polylogarithm iJ(m; x). m is a list representing 
the index vector. We chose the "m" -notation as the standard notation. It 
is possible to give as argument a vector in the "a" -notation, or even a mix 
between the two notations, the result will be automatically converted to 
the "m" -notation. 

« HPL" 

HPL[{0, 0, 1, 1, 0, -1}, X] 

HPL[{3, 1, -2), x] 
HPL[{10, 3, 0, 4), X] 

HPL[{10, 3, 5}, x] 



• HPLMtoA[m_List] ,HPLAtoM[a_List] convert vectors of the "m"- to the 
"a"- and from "a"- to "m" -notation respectively. Both can convert vec- 
tors which mix the two notations. 



HPLMtoA[ {2, 


-3, 0)] 


{0, 1, 0, 


-1, 0} 


HPLAtoM[{l, 


0, 0, 1, 0, -1, 0)] 


{1, 3, -2, 


3} 


HPLMtoA[{0, 


1, -3, 0)] 


{0, 1, 0, 


-1, 0} 


HPLAtoM[{l, 


3, 0, -1, 0)] 


{1, 3, -2, 


3} 



• HPLLogExtract extracts the logarithmic divergences of the HPL in its ar- 
gument at and 1. The result is displayed as function of log(a;), log(l — x) 
or H{l;x),H{0;x) depending on the option settings (see section [Tl.3|) . 



$HPLAutoConvertToKnownFunctions = 


False; 




HPLLogExtract [HPL [{1, 2), x] ] 






HPL[{11, x] HPL[{2}, x] - 2 HPL[{ 


2, 1), X] 




$HPLAutoConvertToKnownFunctions = 


True; 




HPLLogExtract [HPL[{1, 2), x] ] 






-Log[l - x] PolyLog[2, x] - 






2 (| Log[l - x] ^ Log[x] + 






Log[l - x] PolyLog[2, 1 - x] 


- PolyLog[3, 1 


- X] + 


Zeta[3] ) 







• HPLConvertToKnownFunctions returns its argument with HPL's replaced 
by their representation in terms of more common functions, whenever pos- 
sible. This is only needed if the option $HPLAutoConvertToKnownFunctions 
is set False (see section Hi .3(1 . 

$HPLAutoConvertToKnownFunctions = False; 
HPL[{2, 2), X] 

HPL[{2, 2), x] 

HPLConvertToKnownFunctions [%] 

I PolyLog[2, x]^ - 2 PolyLog[2, 2, x] 



• HPLProductExpand returns its argument where the products of HPL's of 
weight wi, . . . ,Wk are replaced by their representation as a linear combi- 
nation of HPL's of weight wi + ■ ■ ■ + Wk- In order to expand all products, 
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HPLProductExpEoid expands the argument (using Expand), so that terms 
of the form 

H{...;x){H{...;x)+H{...;x) + ...). 



are also replaced. 



HPLProductExpand [HPL [ { 4 , 


1}, t] HPL[{3, 2), t]] 




HPL[{3, 2, 4, 1), t] +2 


HPL[{3, 3, 3, 1}, t] + 




HPL[{3, 4, 1, 2), t] + 


5 HPL[{3, 4, 2, 1}, t] 


+ 


12 HPL[{3, 5, 1, 1), t 


+ HPL[{4, 1, 3, 2), t] 


+ 


HPL[{4, 2, 2, 2), t] + 


3 HPL[{4, 2, 3, 1), t] 


+ 


5 HPL[{4, 3, 1, 2), t] 


+ 14 HPL[{4, 3, 2, 1), 


t] + 


27 HPL[{4, 4, 1, 1), t 


+ 4 HPL[{5, 1, 2, 2), 


t] + 


14 HPL [{5, 2, 1, 2), t 


+ 26 HPL[{5, 2, 2, 1}, 


t] + 


36 HPL[{5, 3, 1, 1}, t 


+ 30 HPL [ { 6, 1, 1, 2 1, 


t] + 


40 HPL [ { 6, 1, 2, 1 } , t 


+ 30 HPL[{6, 2, 1, 11, 


t] 



HPLConvertToSimplerArgument returns its argument with HPL's of ar- 
guments ~x, x^ , 1 — xjix — 1) and (1 — a;)/(l + x) replaced by their 
expansion as a sum of HPL's of argument a;, as described in section|Hl 

HPLConvertToSimplerArgument [HPL [ { 2 , 1 ) , -x] ] 

HPL[{-2, -1), x] 

HPLConvertToSimplerArgument [HPL [{2, 1), 1-x]] 

-HPL[{1, 0, 0), x] + Zeta[3] 

HPLConvertToSimplerArgument [HPL [ { 2 , 1 ) , ^ ~ ^ ] ] 

HPL[{~1, -2), x] + HPL[{1, -2}, x] - 

HPL[{-1, -1, -1}, x] + HPL[{-1, -1, 0), x] - 
HPL[{-1, 0, 0), x] - HPL[{1, -1, -1), x] + 
HPL[{1, -1, 0), x] - HPL[{1, 0, 0), x] + 
HPL[{-1, -1), x] Log[2] - HPL[{-1, 01, x] Log[2] + 
HPL[{1, -11, x] Log[2] - HPL[{1, 01, x] Log[2] - 

I HPL[{-1), x] Log[2]^ - i HPL[{1), x] Log[2]^ + 

Zeta [3] 



HPLReduceToMinimalSet returns its argument with the HPL's projected 
to the minimal set, as described in section [S] 

HPLReduceToMinimalSet [HPL[{1, 2, 3}, x] ] 

HPL [{2}, x]^ - 4 HPL [{3), x] HPL [{2, 1), x] - 

HPL [{2), x] HPL [{2), x]2 - 2 HPL [{3, 1), x] ] + 

HPL[{1), x] (HPL[{2), x] HPL[{3), x] - 

3HPL[{3, 2), x] -6HPL[{4, 1), x] ) + 
2HPL[{2, 1, 3), x] -2HPL[{2, 3, 1), x] 



HPLAnalyticContinuation returns its argument with HPL's replaced by 
their analytic continuation. The arguments of the HPL's are taken to be- 
long to the interval specified by the option HPLAnalyticContinuationRegi' 
which can be either 

minftoml the interval — oo to —1 
mltoO the interval —1 to 
onetoinf the interval 1 to cx). 

The HPL's are replaced by their representation in terms of HPL's of ar- 
gument in the interval 0, 1. The choice of the side of the branch cut from 
which the argument is approached is set by the option 

HPLAnaly t i cCont inuat i onS ign 

which can take values —1, 1 or any symbol. 
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HPLAnalyticContinuation [HPL [ { } , 


yl- 




AnalyticContinuationRegion -> mltoO , 




AnalyticContinuationSign -> 1] 






I Pi + HPL [ { } , -y] 






SimpXify [HPXiAnaJ.y'CicCon'tinua'tion 




y] ' 


AnalyticContinuationRegion -> 


onetoinf , 




Analyt ioOont inuat ionS ign 


111 




. 2 

\ iTrTrni ^ riri 




+ 


^-^2 nr^ J_i [ X ± ) , ^ ^ ^ J nt- J_i 


1 + y ^ 




HPL[{1, 1}, y] - I Pi Log 

1 + y 


[2] - 




HPL[{0), yVV^I Log[2] - HPL[ 




Log [2] - 


Log [2] ^ 

2 







If the option HPLAnalyticContinuationRegion is omitted, and if the ar- 
gument is numerical, HPLAnalyticContinuation wiU use automaticaUy 
the appropriate setting. If the option HPLAnalyticContinuationSign is 
omitted, HPLAnalyticContinuation wiU use the value stored in the vari- 
able $HPLAnalyticContinuationSign which is set by default to 1^. 

HPLAnalyticContinuation [HPL [ { 1 ) , |] , 
AnalyticContinuationSign -> -1] 

-I Pi - HPL[{0}, i] - HPL[{1), |] - Log[2] 

Simplify [HPIAnalyticContinuation [HPL [{1, 0}, j^] ] ] 

, 2 

~- + HPL[{1), _ HPL[{1, -1), |] + HPL[{1, 1), |] + 

I Pi Log[2] + HPL[{-1), |] Log[4] + 
HPL[{1), |] (I PI + 2 HPL[{-1), |] + Log[4]) 

It is to be noted that the Mathematica conventions for the analytic con- 
tinuation are not always the same as that of the HPL package. This is 
illustrated by the following example 



HPL[(1}, 


1 


9] 


0.105361 


+ 


3.14159 I 


-Log[l - 


1 


9] 


.105361 




3.14159 I 



Since the substitution of the HPL's through more common functions has 
precedence over the analytic continuation, the option 

$HPLAutoConvertToKnownFunctions 

can interfere with the analytic continuation, as shown in the following 
example. 

$HPLAutoConvertToKnownFunctions = True; 

HPLAnalyticContinuation [HPL [ { } , t ] , 
AnalyticContinuationSign -> 5, 
AnalyticContinuationRegion -> mltoO] 

Log[t] 

$HPLAutoConvertToKnownFunctions = False; 

HPLConvertToKnownFunctions [HPLAnalyticContinuation [ 
HPL[{0}, t] , AnalyticContinuationSign -> 5, 
AnalyticContinuationRegion -> mltoO]] 

I PI 5 + Log[-t] 

This example shows that with the option 
$HPLAutoConvertToKnownFunctions 

set True we loose the control over the sign of the imaginary part (as in 
the first case it will now depend on the Mathematica conventions). 

• MZV [m] is the Multiple Zeta Value (see section llOf) corresponding to the 
index vector m. Their value in terms of mathematical constants are tab- 

^This is the same convention as [41], but opposite to that of [43] 
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ulated® up to weight 8 and systematically replaced. For higher weights, 
the cases covered by Appendix 1X1 are also replaced. 



Table [MZV[ {n} ] , {n. 


2, 8}] 






^ , 2 .4 

Pi Pi 




945' ^e^^L/J, 


9450 ' 


Zeta [ 5 ] , 


MZV[{-2, -2}] 
















MZV[{3, 1, 3, 1, 3, 


1)] 






„ . 12 
Pi 








43589145600 









• The function $HPLOptions gives a list of the options of the package and 
their current values. 

$HPLOptions 

{ SHPLAnalyticContinuationSign -> -1, 

$HPLAutoConvertToKnownFunctions True, 
$HPLAutoConvertToSlmplerArgument -> False, 
$HPLAutoLogExtract -> False, 
$HPLAutoProductExpand -> True, 
$HPLAutoReduceToMinimalSet -> False} 



11.2 Functions modified 



• We defined the derivatives of the HPL's as described in section |3 The 
integration showing up in the recursive definition of the HPL's is also 
implemented. 



$HPLAutoConvertToKnownFunctions = False; 
D[HPL[{2, 10, 4), X], X] 
D[HPL[{-1, 4, -6}, X], X] 
D[HPL[{1, 3, 7}, X], X] 

HPL[(1, 10, 4), x] 

X 

HPL[(4, -6}, x] 



HPL[(3, 7}, x] 

1 



Integrate [ 
Integrate [ 



HPL [ { 2 , 


4), 


X] 


1 - 


X 




HPL [ { 2 , 


4), 


X] 


X 






HPL [ { 2 , 


4), 


X] 


1 + 


X 





Integrate [ 

HPL[{1, 2, 4}, t] 
HPL[{3, 4), t] 
HPL[{-1, 2, 4), t] 



{X, 0, t)] 
{X, 0, t}] 
{X, 0, t}] 



The function Series is able to expand HPL's around x = and x = 1. 



Series [HPL [{2, 3, 4, 5}, x] , {x, 0, 7)] 

x^ 193 x^ 1026311 x^ 241822151 x^ 

6912 * 1382400 * 8957952000 * 2633637888000 ^ 
SeriesData [Global 'x, 0, {), 4, 8, 1] 

Series [HPL [{1, 2, 3), x] , {x, 1, 1}] 

SeriesData[l - Global "x, 0, 

, -143 Pi^ + Pi^ I.og[l - Global'x] Zeta[3] y^r^f^i' 
* 45360 3 ZetaL3J 



9 Log[l 
. 4 



Global'x] Zeta[5] 

2 



-Pi 
72 



2 Zeta[3] 



Log[l - Global'x] Zeta[3]), 



^these tables are those of the FORM package harmpol [39] 
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11.3 Working with options 

The package HPL has some options to control its behavior. They set the preferred 
form in which expressions are displayed. The option can be locally overridden 
by the functions described above. The effects of the options are described in 
the following. 

$HPLAutoLogExtract: If True the logarithmic divergences log(l — a;) and log(a;) 
are automatically extracted from the HPL's following the procedure de- 
scribed in section[S| These divergences are displayed as log or as H{Q/1\ a;), 
depending on the value of the option $HPLAutoConvertToKnownFunctions. 
The default setting is False. 

$HPLAutoProductExpand: If True the products of HPL's arc automatically con- 
verted into a sum of HPL of weight equal to the sum of the weights of the 
two factors, as described in section^] Default is False. 

$HPLAutoConvertToKnownFunctions: If set True, the HPL's will be converted 
to more common functions (logs, poly logarithms, polynoms, Nielsen poly- 
logs) if possible, using the identities of Appendix lUl This might be coun- 
terproductive when the properties of the HPL's are more explicit in the 
HPL form than in their equivalent representation. Furthermore, if this op- 
tions is set True while using the analytic continuation described above, the 
result may be wrong, as Mathematica does not have different conventions 
for the analytic continuation. Default is False. 

$HPLAutoReduceToMinimalSet: If set True, the HPL's will be automatically 
reduced to a minimal basis (up to weight 8). This only makes sense if one 
does not expand the obtained products again, or if the factors of smaller 
weight can be replaced by their expression in terms of known functions. 
Therefore, for the reduction to be performed, one has to have the option 

$HPLAutoProductExpand equal False or 

$HPLAutoConvertToKnownFunctions equal True. 

If this is not fulfilled, the option will have no effect. Default is False. 

$HPLAutoConvertToSimplerArgument: If set True, the HPL's of arguments 
—X, , 1 ~ X, l/x, x/{x — 1) and (1 — x)/{l + x) will be automatically 
substituted by their representation in term of HPL's of argument x along 
the lines of section|Hl Default is False. 

$HPLAnalyticContinuationSign: If set to 1 the analytic continuation of the 
HPL's is taken for arguments with an positive infinitesimal imaginary 
part, if set to -1 with an negative one. This is only the default setting and 
can be overridden by specifying the option AnalyticContinuationSign 
in the function HPLAnalyticContinuation, as described above. Default 
setting is -1-1. 

11.4 Numerical evaluation 

If the argument of the HPL is numerical, the package evaluates its value and 
gives the result in the same precision as the precision of x. In [1], the authors 
propose to use the argument transformation x {l—t)/{l+t) for the evaluation 
of HPL's with argument near one or negative with a large absolute value. This is 
only of advantage if one wants to evaluate the whole set of HPL's for this given 
X. Since for low weight, we have a representation in terms of usual functions 
(which Mathematica can evaluate fast and precisely), this change of variable is 
not useful in this case. For high weights, the expression of the HPL's in term 
of HPL's of argument is quite large, so that the time gain in the conver- 
gence is compensated by the number of different HPL's to evaluate. Therefore, 
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we used the series expansion systematically for the evaluation of the numerical 
value of the HPL. The convergence is not very good for values near 1, where not 
all digits displayed are accurate. There is a FORTRAN package for the precise 
numerical evaluation of HPL's [41]. 



HPL[{2, 1, 


4, 


5), 


1] 




MZV[{2, 1, 


4, 


5)] 






HPL[{2, 1, 


4, 


5), 


1. 


] 


. 0272374 










HPL[{2, 1, 


4, 


5), 






HPL[{2, 1, 


4, 


5), 


b 




HPL[{2, 1, 


4, 


5), 


2. 


00000000000000000000000000000^ 


.00017490 


728 


2457760 


677740168010020 


N[HPL[{2, 


1, 


4, 5), 


|], 50] 


.00017490 


728 


2457760 


67 774016801007 622686190198808601431 



For argument outside the interval 0, 1, one can give the sign of the infinitesimal 
imaginary part to be taken in account for the analytic continuation, as for the 
function HPLAnalyt icCont inuat ion 

HPL [{2, 1, 0), 1.4, AnalyticContinuationSign -> 1] 

-1.99574 - 1.73879 I 

HPL [{2, 1, 0), 1.4, AnalyticContinuationSign -> -1] 

-1.99574 + 1.73879 I 

For the numerical values of the MZV, we implemented the procedure described 
in [43,44]. The numerical values of the constants se, sra, sn, ssa, ssh, sgc and ssd 
correspond to the results given by the EZf ace application [47]. 

N[MZV[{4, 3)]] 

.0851598 

N[MZV[{2, 5, 4, 3)], 40] 

0.0000 97 163 968999849198015543839 97 828478502511 

Due to the possibility of evaluating the HPL numerically. Plot is able to repre- 
sent HPL's 



Plot [{Re [HPL [{- 


-1, - 


-1, - 


1, 1), t]], 




Im[HPL[{-l, 


-1, 
Pi 


-1, 


^'^}, (t, -5, 


5), 


PlotStyle -> 


{Hue[l], 


Dashing [{0.02, 


02)] ), 


PlotRange -> 


{-4 

4 
3 


4)] 






-4- ^ -2 \ 


-1 




2 4 






-2 










-3 








- Graphics - 


-4 









We checked numerical agreement with [41] at double precision accuracy. 

12 Conclusion 

In this paper, we presented an implementation of the harmonic polylogarithm of 
Remiddi and Vermaseren [1] for Mathematica. It contains an implementation of 
the product algebra, the derivative properties, series expansion and numerical 
evaluation. The analytic continuation has been treated carefully, allowing the 
user the full control over the definition of the sign of the imaginary parts. Many 
options enable the user to adapt the behavior of the package to his specific 
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problem. This package is already used by the package HypExp [36] to treat 
the HPL's and MZV's appearing in the expansion of hypergeometric functions 
around integer-valued parameters. 
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A Multiple Zeta Values 

We list here some of the identities found in [44, 45] . 

C(2,l) = C(3) (43) 

2C(m,l) = mC{m + l)-^C{m~k)C{k + l), 2 < m & Z (45) 

fe=l 

C(2,{in = an + 2) (46) 
C(3,{in = C(^ + 2,l) (47) 

-J- 9 1 " 

= ^-((n + 3)--^£^('' + mn + 2-k) (48) 



(2n-t- 1)! V2 

= m^aT" (50, 



(4n-|-2)! V2 

2-77^" 

C({3,in = 4-"C({4r) = ^^^^ (53) 

n 

C(2,{l,3r) = 4-"^(-l)'=C({4r){(4A; + l)C(4fc + 2) 

k 

-4^C(4i-l)C(4fc-4i3)} (54) 
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A.l Colored MZV 

Here we list some identities for MZV with negative arguments found in [45] 

. ((1 + V2 +(1-V2 )))(56) 

C({-6}„) = (^;7^2(^ + ^'"^'("^)"^"'''^^' 

^\ 6n+3 / ^\ 6)1+3 

«{-!)..) = (-irsn^ 

where the sum in the last equation is over all non negative integers satisfying 
Sfe>o = n 

C(-1,{1}„) = (59) 

n! 

C(-l,-l,{l}n) = -L*„+2(l/2) (60) 

B Conversion of HPL at unity to MZV 

We first prove 

H{-m;x) = -C{-m)=Y^ — — — 

i=l 



oo oc 

OO ^ oc - 



4=1 



= -2-™C(m) + (C(m) - 2-'"C(™)) 

= (1 - 2i-'")C(m) (61) 

For the proof of the formula for higher depth we start with the definition p2|l 
of the coefficients of the series expansion of the HPL and first prove that 

(-l)'+iZ,(mi,...,fc) - (-l)^+iZ,(-mi,...,fc). (62) 
We show that it holds for depth 1, 

(-i)-*z.(™,) = (-i)-^^^^^^ 

i^iy+^-^ = Z,{-m,), mi>0 

(63) 

{-iy+^ ^ ' =-. =Z,(-TOi), mi<0. 
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For mi,...,fe with mi > we get 

{-iy+'Z,{mi_k) = 



12 = 1 



^(_l)i.+i(_l)^Z,,(-m2,...,fe) 



i2=l 

(-l)'=+iZi(-mi,...,fc), 



(64) 



and for mi < 0, 



(-l)'+iZi(mi,..., 



(1 



^|mi| 



i-1 



5^(-ir+iZ,,(m2,...,fe) 



i2 = l 



j|mi| 



5^(-l)^Z,,(-m2,...,fc) 



i2 = l 



= (-l)'=+^Z,(-mi,...,fc). 
The next step is to prove for A; > 1 

sgn(TOi)' %gn(mim2)*^ sgn(mfe_imfc)*'= 



(65) 



12=1 ik 



(66) 



where N{m\^,,,^k) is given by (—!)#" with #n the number of negative indices 
in the index vector m. 

The proof is again through induction in the depth. We first test the claim 
for depth 2. For mi positive 

1 ^ „ . . 1 ^ sgn(m2)*=+i 



^i(mi,m2) = ^ 5^ ^^2(^2) = 7^ E 



i2 = l 



i2=l 



j|m2| 



1 '"^ 1 

— E — 

jmi Z_-^ jm2 

12 = 1 



1 ^ (-1)^^ + 1 

12 = 1 



1 

= -— E 



(-1)^ 



^mi ^|m2| 

12 = 1 



m2 > 



m2 < 



(67) 



which is the expected result. For mi negative one has 

(-1)' ^ . ^ (-1)^ sgn(m2)*^+i 



Zj(mi,m2) 



E^^^(-^) = SE 



i2 = l 



i2 = l 



i\m2\ 



{-ly 



E — 



22 = 1 



(-1)' (-1)'^ + ! 

j|mi| l\m2\ 

12 = 1 



hi): 

j\m2\ 



12=1 



m2 < 



m2 > 



(68) 
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Wc turn now to the general case k > 2. For mi positive, one has 
1 '"^ 

^i(mi,...,fe) = — X! ^»2("i2,...,fc) 

^ Z— 1 ik-± 



^ ' sgn(m2)'^ sgn(TO2™3)''' sgn(mfc_imfc)"'= 

^mi ^""^^■■-'^^ Zl^ Zl^ -Imal -Imgl .|mfc| 

12 = 1 ifc *2 ^t- 



and for mi < 



_ sgn(mi)' !^ 

12 = 1 

42 = 1 Ifc 

-sgn(m2)'" sgn((-m2)(-TO3))'" sgn((-mfc-i)(-mfc))"' 

• Imal •|m2| '■' -Irrifcl 

'2 '2 *fc 

22 = 1 ifc *2 *fc 



since 

(-l)'=7V(-m2,...,fc) = 

— (^_2^')-#n("i2,...,fc-l) — (■_2^')#n(™l,...,fc) ^72) 

The connection between the HPL and the MZV is easy to figure out from the 
definition ^ of the MZV. One obtains the result ^ 

H{{mi,...,k}, 1) = iV(mi,...,fc)C(mi.....fe). (72) 

C Table of representation through more com- 
mon functions 

C.l Weight 2 

H{{l,l};x) = ^log(l-a;)' (73) 
if ({1,-1}; a;) = Li2((l - a;)/2) - log(2) log(l - a;) - iz2 (74) 

i7({-l,l};a;) = Lz2((l + 2:)/2) - log(2) log(l + x) - iz2 (75) 

if ({-1,0}; a;) = log(l + a:) log(x) + Li2(-x) (76) 
H{{-2};x) = -Li2{-x) (77) 
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C.2 Weight 3 



H{{l,2};x) = -2Si,2{x)~\og{l-x)Li2{x) (78) 

H{{l,l,l};x) = -ilog(l-a;)' (79) 

H{{2,l};x) = Sr,2{x) (80) 

^.({2,-1};.) = Li.i^J-U,[^)-U. ^'^^ 

-Li^ix) + log(l + x)Li2{^) + log(l + x)Li2{x) 

+ilog(2)log(l+a;)2+Lz3(^i) (81) 

H{{-2,l}-x) = -S,^,{x)+Lis(^^^-Us(^^^-U-,{-x) 

+Lz3 + Li3{x) + log(l - x)Li2{-x) 

+ log(l - x)Li2 - log(l - x)Li2{x) 

+ i log(2) log(l - xf - l/61og(l - xf (82) 
if ({-2,-1}; a;) = Si,2{-x) (83) 
({1,-1,-1}; a.) = -Uog{^\\og{l + xf-uj\ 



-\og{l+x)Li2(^-±^^+Li,(^-±^ 



(84) 



i?({-l,l,l};x) = -log(^^jlog(l-x)^ + L^3(^2j ^^^^ 

+ log(l - a;)Li2 (^^) - L,;3 (^^) (86) 

if ({-1,-2}; a;) = H {{-2}-x) H {{-l}-x) - 2H {{-2,-l}-x) (87) 
if({-l,2};x) = H{{-l]-x)H{{2}-x)-H{{-2,l}-x) 

-H{{2,~l}-x) (88) 
if({-l,l,-l};x) = iJ({-l};x)i?({l,-l};x)-2i?({l,-l,-l};a;) (89) 

C.3 Weight 4 

if ({1,3}; a;) = ~Li2{xf -\og{l - x) Lh{x) (90) 

//({l,l,2};x) = ^\og{l-xf Li2{x)+2\og{l-x)Si,2{x) 

+35i,3 (x) (91) 

F ({1,2,1}; a;) = - log (1 - a;) 5i,2 (a;) - 35i,3 (x) (92) 

if ({2, 2}; a;) = -282,2 {x) + ^Li2{xf (93) 

if ({3,1}; a;) = ^2,2 (x) (94) 

if ({2,1,1}; x) = 5i,3(a;) (95) 

if ({1,1, 1,1}; a;) = -Llog(l-ar)^ (96) 
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C.4 Arbitrary weight 



H{{n};x 
H{{n};x 

//({"!}; X 

i7(r(-l)};x 
H{{nr'l};x 



Lin (x) , n > 

/ ^Nn log(l-a:)" 
n! 

log(l + a:)" 
n! 

'5'n— l,p (3^) 



(97) 
(98) 

(99) 

(100) 
(101) 
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